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Q ] Abstract 

I As shown in Ref. [11], the Fermi nuclear matrix element Mp^ of neutrinoless double beta 

decay can be reconstructed if one is able to measure the isospin-forbidden single Fermi transition 
matrix element from the ground state of the final nucleus to the isobaric analog state (IAS) of 
the initial nucleus, for instance by means of charge-exchange reactions of the (n,p)-type. Here, 
O ■ simple estimates for ^^Se are made which show that indeed the tiny cross section anp{0~j — t- IAS) 

^ . is dominated by the admixture of the double IAS in the ground state of the final nucleus provided 

that the isospin mixing is weak and can be treated perturbatively. A measurement of such a cross 
section would definitely be a very difficult task, but it can advance a lot our knowledge of the Oz//3/3 
^ ■ nuclear matrix element. 

^ I Neutrino is the only known spin-| fermion which may be truly neutral, i.e., identical with its own 
] antiparticle. In such a case one speaks about Majorana neutrino, to be contrasted with Dirac neutrino 
CN which is different from its antiparticle pQ. Majorana neutrinos naturally appear in many extensions of 
^ the standard model (see, e.g., p]). 

T— I . A study of nuclear neutrinoless double beta (Oz//3/3) decay z^\n — > z+2^^n-2 + 2e~ offers the 
j> ! only feasible way to test the charge-conjugation property of neutrinos. The existence of 0^(3(3 decay 
, would immediately prove neutrino to be a Majorana particle [3l H]. The decay also allows to probe the 
^ ■ absolute neutrino masses at the level of tens of meV. 
■ — ■ The next generation of 0i//9/9-decay experiments (see, e.g., Ref. ^ for a recent review) has a great 
discovery potential. Accurate knowledge of the relevant nuclear matrix elements M'^^ will be crucial to 
reliably deduce the effective Majorana mass from the future experimental data. 

Several theoretical approaches have been used to evaluate M^" . The present world status of the 
results on M^'^ for the light neutrino mass mechanism is shown in Fig. [1] (adopted from Ref. [5]). Here, 
we are not going into a detailed discussion of advantages and disadvantages of different models which 
can be found in Ref. 0. One only notices that the calculated NSM M^" scatter substantially, up to a 
factor of 2. 

In such a situation, it would be extremely important to find a possibility to determine M^'^ exper- 
imentally. There have been attempts to reconstruct the nuclear amplitude of two-neutrino decay 
from partial one-leg transition amplitudes to the intermediate l"*" states measured in charge-exchange 
reactions [10]. However, such a procedure can consistently determine M^^ only if a transition via a 
single intermediate l"*" state dominates M"^", since relative phases of different contributions cannot be 
measured. Trying the same way to reconstruct M"^ seems even more hopeless, since many intermediate 
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Figure 1: (Color online) Neutrinoless double beta decay transition matrix elements M^*^ cal- 
culated in different approaches: the quasiparticle random phase approximation (QRPA) 
the nuclear shell model (SM) [7], the projected Hartree-Fock-Bogoliubov method (PHFB) [8J 
and the Interacting Boson Model (IBM-2) [9J. The error bars in the QRPA results are calcu- 
lated from the highest and the lowest values of M^'^ obtained in the calculations with three 
different single-particle basis sets, two different axial vector coupling constants qa = 1-25 
and gA = 1-00 (the quenched value), and two different treatments of short range correlations 
(Jastrow-like and the UCOM). 



states of different multipolarities (with a rather moderate contribution of the 1^ states) are virtually 
populated in the Oi//3/3 decay due to a large momentum of the virtual neutrino. 

An alternative proposal was put forward in a recent work [TT]. It exploits the similarity between the 
Fermi part of the neutrino potential in Oz//3/3 decay and the radial dependence of the two-body Coulomb 
interaction. The latter is well-known to be the leading source of the isospin breaking in nuclei [12]. As 
shown in Ref . [11] , the Fermi part M^^ of the total matrix element M^'^ can be related to the Coulomb 
mixing matrix element between the ideal double isobaric analog state (DIAS) of the ground state (g.s.) 
of the initial nucleus and the g.s. of the final nucleus. This ideal DIAS would be an exact nuclear state 
if the isospin symmetry were exact. As a result of the Coulomb mixing, the single Fermi transition 
matrix element {Of\T~\IAS) between the isobaric analog state (IAS) of the g.s. of the initial nucleus 
and the g.s. of the final nucleus becomes non-zero as is illustrated in Fig. [2l Thus, having measured 
this single Fermi transition matrix element, e.g., by charge-exchange reactions, the Oz//3/3-decay matrix 
element Mp" can be reconstructed. 

The master relation, derived in Ref. [H] by making use of the isospin symmetry of strong interaction, 
represents the matrix element Mp^ in the form of an energy-weighted double Fermi transition matrix 
element: 

M'/ = -^,T.^s{0f\f-\0t){0t\f~m. (1) 

° s 

Here, = J2a '''a is the isospin lowering operator, the sum runs over all O"*" states of the intermediate 
{N — 1, Z + 1) isobaric nucleus, ujs = Eg — {Eq. + Eq^)/2 is the excitation energy of the s'th intermediate 



2 



IT T -2> 
^IDIAS> 




Figure 2: (Color online) Double Fermi transition from the ground state |0j) of the initial 
nucleus would go exclusively to its DIAS in the final nucleus if the isospin symmetry were 
exact. The isotensor part of the Coulomb interaction admixes this ideal DIAS to the g.s. 
|0/) of the final nucleus. Because of the mixing, the single Fermi transition matrix ele- 
ment {IAS\T~^\Of) becomes non-vanishing, and could be measured, e.g., by charge-exchange 
reactions. 



state relative to the mean energy of g.s. of the initial and final nucleus. 

As argued in Ref . [H] , the expression ([T]) must be dominated by the amplitude of the double Fermi 
transition via the IAS of the initial g.s. into the final g.s.: 

M^-^ ^ -^ujAs{Of\f-\IAS){IAS\f-\0,). (2) 

Here, the second Fermi transition amplitude is non-vanishing due to an admixture of the ideal DIAS 
\DIAS) = , |0+) in the g.s. of the final nucleus: {Of\f-\IAS) = {Of\DIAS){DIAS\f-\IAS), 

y 4To(2To-l) 

To = {N — Z)/2 is isospin of the g.s. of the initial nucleus. 

Therefore, the total M'^" can be reconstructed according to Eq. ([2]), if one is able to measure 
the AT = 2 isospin-forbidden matrix element {IAS\T~^\Of) , for instance in charge-exchange reac- 
tions of the (n,p)-type. Using recent QRPA calculation results for Mp" , this matrix element can 
roughly be estimated as {IAS\T^\Of) ~ 0.005, i.e. about thousand times smaller than the first-leg m.e. 
{IAS\T~\Oi) ~ y/N — Z. This strong suppression of {IAS\T^\Of) reflects the smallness of the isospin 
mixing effects in nuclei. 

The IAS has been observed as a prominent and extremely narrow resonance, and its various features 
have well been studied by means of {p,n), {^B.e,t) and other charge-exchange reactions on the mother 
nucleus. In this case the reaction cross section at the zero scattering angle can be shown to be propor- 
tional to a large Fermi matrix element {IAS\T~\Oi) ^ ■\/27o [IS]. Extraction of a strongly suppressed 
matrix element {IAS\T'^\Of) from the tiny cross section of the (n,j9)-type charge-exchange reactions 
on the final nucleus might only be possible if there is a proportionality similar to the {p, n) channel. 
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Therefore, a realistic analysis of the corresponding reaction mechanism is needed. 

First of all, it is easy to argue that isospin of the projectile may not be larger than T = 1/2. Indeed, 
the main components of the wave functions |0/) and 1/^5") have total isospin different by two units. 
Therefore, already for a projectile with isospin T = 1 there exists a common entrance and exit isospin 
channel arising from recoupling of the isospin T = 1 of the projectile with the main components of the 
wave functions of the target and daughter nuclei. In such a case extraction of the information about 
small isospin impurities from the corresponding reaction cross section seems barely possible. 

Thus, the only feasible probes must be of isospin T = 1/2 {{n,p), (t,^He),. . .reactions). However, 
even for these probes it is still not guaranteed that the reaction cross section is proportional to a strongly 
suppressed matrix element {IAS\T~^\Of) since the other isospin impurities in the wave functions |0/) 
and \IAS) may have larger effect on the reaction cross section. In the following we shall assess the case 
of the {n,p) reaction at the zero scattering angle. We also restrict ourselves to the case of a single, 
well-isolated, IAS as it appears in rather light nuclei like ^^Ca, neglecting therefore the spreading of 
the IAS that becomes important in heavier nuclei. In such a case the Coulomb mixing can be treated 
perturbatively that significantly simplifies the consideration. 

The g.s. of the initial nucleus has the largest isospin projection = Tq among all the nuclei 
involved, therefore it can be considered as a state of the pure total isospin T = Tq: \0f) = |ToTo)i. 
However, small isospin admixtures induced by the Coulomb mixing should be taken into account in the 
wave functions of the IAS and the g.s. of the final nucleus: 

\IAS) = \ToTo-l)^ + J2(^s\To-lTo-l)s, (3) 

s 

\0p = \To-2To-2)f + J2Ps\To-lTo-2)s + ^\ToTo-2)i. (4) 

s 

It suffices for description of the IAS wave function Eq. to take into account only the admixtures 
of all 0"*" states with the normal isospin Tq — 1 in the intermediate nucleus. For the g.s. of the final 
nucleus, Eq. (jl]), the admixtures of all T-greater 0+ states must be considered as well as the only 
AT = 2 admixture of the ideal DIAS |ToTo — 2)j whose amplitude 7 is the quantity to be determined 
experimentally. The amplitudes a and P of the AT = 1 isospin admixtures in Eq. fl3ll^ are mainly 
determined by the mean Coulomb field Uc, while the amplitude 7 is essentially determined by the 
isotensor two-body Coulomb interaction. 

In order to be able to determine 7 experimentally, one would need the reaction channel going through 
the admixture of the ideal DIAS: 7 |n) ® |To Tq - 2)^ ^ \p) ® [Tq Tq - 1)^, with T = Tq ± i, T^ = Tq - |, 
to dominate the cross section of the {n,p) reaction at the zero scattering angle. In fact, the AT = 1 
isospin admixtures allow for other open competitive channels: |n) ® /3|To — 1 Tq — 2) — )■ \p) ® |To Tq — 1) 
(with T = To — |, channel I) and \n) (g) |To — 2 Tq — 2) — )• \p) eg) «|To — 1 Tq — 1) (with T = Tq — |, channel 
II), whose relative contributions we are going to estimate. 

Let us first estimate the cross section in the channel I. The corresponding T-matrix for the one-step 
reaction is: 

Tm(0+ ^ IAS) oc {p\ 0, (To To - l|Kt.|n) ® J] /3.|To - 1 Tq - 2)„ (5) 

s 

where T-greater states in the final nucleus are the analog states of 0"'' states in the intermediate nucleus: 
|To — 1 To — 2)s = ^2To-2 1-^0 ~ 1 ^0 ~ and Vstr is the interaction between the projectile and the target 
which is assumed to be dominated by strong interaction. The admixture amplitudes Ps in the first order 
of perturbation theory with respect to the Coulomb mean field read: 

^ _ .(To-lTo-2|[/c|To-2To-2); _ ,(To - 1 To - l|?7+|To - 2 Tq - 2)^ 

^2iTo-l){E, + Ac-E,^) 



E. 
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where the charge-changing Coulomb mean field is defined as = [T"*", Uc] = ~ 



R) inside the nucleus, and Ac is the Coulomb displacement energy. Then, by making use of isospin 
symmetry of Vgtr the expression (|5]) can further be transformed to acquire the form: 



Tm(o+^M^)(x5: 



^(To-2To-2|;7g|To-lTo-l), 
V2To{E, + Ac - E^+) 



,{To - ITo - l\an{k')Vstral{k)\ToTo),. (7) 



The matrix element s(To — 1 Tq — l|a„(fc')V^str'^p(fc)|^o ^o)« determines the T- matrix for the direct forward- 
scattering {p,n) reaction on the g.s. of the initial nucleus with the excitation of monopole states in the 
intermediate nucleus, mainly those which form the isovector monopole resonance (IVMR). One would 
need a reaction code to calculate such a T-matrix, but for a rough estimate we assume here that this 
{p, n) reaction cross section is proportional to the respective monopole strength of a O"*" component of 
the IVMR: 

ap„(0+ ^ IVMRs) = ao \{IVMRs\R-\Ot) 

where do is the unit cross section, and i?^ = J2a ^ '^a is the monopole charge-changing operator. With 
this simplification one can get the ratio of two cross sections: 



ap„(0+ ^ IVMR) 
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(9) 



The nominator in this expression has a form of the squared amplitude of a double monopole transition 
between the g.s. of the initial and final nuclei, and the denominator contains the total /3~ monopole 
strength of the IVMR excited from the g.s. of the initial nucleus. 
The same transformation with the T-matrix in the channel II 



Ti? oc, (To - 2 To - 2K{k') aj (fc) \T, -IT,- l|^c-|^oTo). 



leads to the second ratio: 
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(10) 



where the denominator contains the total (5^ monopole strength of the IVMR excited from the g.s. of 
the final nucleus. 

Taking as an estimate anpi'jDIAS — )■ IAS) ~ 10~^(7p„(0j^ — )■ IAS) for the (n,p)-reaction cross 
section at the zero scattering angle which is due to the admixture of the DIAS, one arrives at the 
following ratios of the cross sections in question: 



a(5(0+^M^) 



anpi'yDIAS IAS) 



ag)(0+->M5) ^ ^o6^(ii) anpiOj^IVMR) 



anpi-fDIAS IAS) 



ap^jOt ^ IVMR) 
apniOf ^ IAS) 
r„p(0+ ^ IVMI 
apniOf ^ IAS) 



(12) 
(13) 
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For a numerical estimate we choose ®^Se as a representative medium-mass /3/3-decaying nucleus. Both 
AN = and AA^ = 2 particle-hole (p-h) excitations contribute to the monopole matrix elements (the to- 
tal monopole strength is dominated by the contribution of AA^ = 2 p-h excitations). For a simultaneous 
description of both types of excitations a very large single-particle (s.p.) basis should be used in a nuclear 
model which at the same time should preserve isospin. To circumvent this difficulty, here two separate 
methods are used for the estimate: AA^ = 2 p-h excitations are treated in the independent Bogolyubov 
quasiparticle picture within a large, 7hu!, s.p. space, whereas AA^ = p-h excitations are described 
within the continuum- QRPA with zero-range Landau-Migdal forces. The estimates of the contribution 
of the AA^ = 2 p-h excitations are: r*^^^ = 6.6-10^^; r^^^^ = 1.3-10^^. In the continuum- QRPA calculation 
of the contribution of the AA^ = p-h excitations the following substitution is used: R~ — ?■ R" — aT~ 
(ct = 77T^!5— r^) which allows to subtract from the monopole matrix elements the contribution of a 
spurious admixture of the IAS in monopole states that is due to a incomplete isospin conservation in 
the model. The estimates of the contribution of these AA^ = p-h excitations are: r*-^^ = 6.4 ■ 10~^; 

r(") = 2.2-10-^ Taking into account the ratios ~^ ~ 0.1, ^"p^^^/ ~^ ~ 0.01 [H 



cr. 



pn 



(J W (Ot ^ IAS) a^") (Ot -> IAS) 

one arrives at the conclusion that — —— — — ~ — —rr. tttt ~ 0.1. So, indeed the 

anp{iDIAS IAS) anpiiDIAS IAS) 

cross section cr„p(Oj — )■ IAS) at the zero scattering angle is dominated by the admixture of the ideal 

DIAS in the g.s. of the final nucleus. 

In the above analysis the assumption of a weak mixing between the IAS and the O"*" states of normal 

isospin in the intermediate nucleus has been employed. However, in heavy nuclei the density of O'^ states 

in the vicinity of the IAS becomes so high that the IAS gets spread (the resonant mixing occurs). This 

effect is still to be taken into consideration. Nevertheless, in the case of the IAS of ^^Ca (T = 4, = 3) 

in ^^Sc the weak-mixing assumption seems to be realized. The corresponding estimates will be published 

elsewhere. 

To conclude, it has been shown in Ref. [HJ that the Fermi OuPP nuclear matrix element can be 
reconstructed if one is able to measure the isospin-forbidden matrix element from the g.s. of the final 
nucleus to the isobaric analog state in the intermediate nucleus, for instance by means of charge- 
exchange reactions of the (n,p)-type. Here, simple estimates have shown that indeed the tiny cross 
section cT„p(Oj — )■ IAS) is dominated by the admixture of the DIAS in the g.s. of the final nucleus 
provided the mixing between the IAS and the 0^ states of normal isospin in the intermediate nucleus 
is weak and can be treated perturbatively. A measurement of such a cross section would definitely be 
a very difficult task, but it can advance a lot our knowledge of the Oz//3/3 nuclear matrix elements. 

The author acknowledges support of the DFG within the SFB TR27 "Neutrinos and Beyond" . 



References 

[1] B. Kayser, F. Gibrat-Debu, and F. Perrier, The Physics of Massive Neutrinos (World Scientific, 
Singapore, 1989); R.N. Mohapatra and P.B. Pal, Massive Neutrinos in Physics and Astrophysics 
(World Scientific, Singapore, 1991). 



[2] R. N. Mohapatra et a/.arXiv:hep-ph/0412099 R. N. Mohapatra and A. Y. Smirnov, Ann. Rev. 
Nucl. Part. Scz. 56 (2006) 569. 

[3] F. Boehm and P. Vogel, Physics of Massive Neutrinos, 2nd ed. (Cambridge University Press, 
Cambridge, 1992); A. Faessler and F. Simkovic, J. Phys. G 24 (1998) 2139. 

[4] Frank T. Avignone HI, Steven R. Elliott, and Jonathan Engel, Rev. Mod. Phys. 80 (2008) 481 . 



6 



[5] A. Escuderos, A. Faessler, V. Rodin, F. Simkovic. ' arXiv:1001.3519l [nucl-th]. 

[6] F. Simkovic, A. Faessler, V.A. Rodin, P. Vogel, and J. Engel, Phys. Rev. C 77 (2008) 045503 . 

[7] E. Caurier, F. Nowacki, A. Poves, Eur. Phys. J. A 36 (2008) 195; E. Caurier, J. Menendez, F. 
Nowacki, A. Poves, Phys. Rev. Lett. 100 (2008) 052503 ; J. Menendez, A. Poves, E. Caurier, F. 
Nowacki,iVMc/. Phys. A 818 (2009) 139 . 

[8] K. Chaturvedi, R. Chandra, P. K. Rath, P. K. Raina, J. G. Hirsch,P/i?/s. Rev. C 78 (2008) 054302 ; 
P. K. Rath, R. Chandra, K. Chaturvedi, P. K. Raina, J. G. Hirsch, Phys. Rev. C 80 (2009) 044303 
R. Chandra, K. Chaturvedi, P. K. Rath, P. K. Raina and J. G. Hirsch, Europhys. Lett. 86, 32001 
(2009). 

[9] J. Barea, F. Iachello,P%s. Rev. C 79 (2009) 044301 . 
[10] D. Frekers, Prog. Part. Nucl. Phys. 64 (2010) 281 . 
[11] Vadim Rodin and Amand Faessler, P/i^/s. Rev. C 80 (2009) 041302(R) . 

[12] N. Auerbach, J. Hiifner, A. K. Kerman and C. M. Shakin, Rev. Mod. Phys. 44 (1972) 48 ; N. Auer- 
bach, Phys. Rep. 98 (1983) 273 . 

[13] T.N. Taddeucci, C.A. Goulding, T.A. Carey, R.C. Byrd, CD. Goodman, C. Gaarde, J. Larsen, D. 
Horen, J. Rapaport, E. Sugarbaker, Nucl. Phys. A 469 (1987) 125 . 

[14] R.G.T. Zegers, private communication. 



7 



